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ABSTRACT 

We have derived the absolute proper motion (PM) of the globular cluster M55 using a 
large set of CCD images collected with the du Pont telescope between 1997 and 2008. 
We find {fiaCosS, fis) = (—3.31 ± 0.10, —9.14 ± 0.15) mas/yr relative to background 
galaxies. Membership status was determined for 16945 stars with 14 < F < 21 from 
the central part of the cluster. The PM catalogue includes 52 variables of which 43 are 
probable members of M55. This sample is dominated by pulsating blue straggler stars 
but also includes 5 eclipsing binaries, three of which are main sequence objects. The 
survey also identified several candidate blue, yellow and red straggler stars belonging 
to the cluster. We detected 15 likely members of the Sgr dSph galaxy located behind 
M55. The average PM for these stars was measured to be {^aCosS, /i5)=(— 2.23±0.14, 
-1.83±0.24) mas/yr. 

Key words: astrometry - globular clusters: individual: NGC 6809 (M55) - binaries: 
eclipsing - blue stragglers - galaxies: individual: Sagitarrius dSph 
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1 INTRODUCTION 

M55 (NGC 6809) is a metal-poor globular cluster (GC) in 
the Galactic halo (1 = 9°, b = —23°), discovered by Nicholas 
Louis de Lacaille in 1752. Over 70 variable sta r s are known 
in the field of the cluster Jciement et al" 19991 : 1 lech et al.l 
ll999l : |Pvch et al.ll200ll : lKaluznv et al.ll2010l . hereafter JKIO). 
About half of these are located in the blue straggler region 
of the color-magnitude diagram (CMD) and 13 are RR Lyr 
type p ulsators. The cluster is close to us {{m~ M)v = 13.87; 
lHarrisllT996i ) and is relatively sparse. This makes it an ideal 
target for a measurement of the cluster proper motion (PM). 
Several authors have shown that for nearby GCs ground 
based CCD observations with temporal baselines of 3-5 years 
are sufficien t to efficiently separate member stars from field 
inter l opers (I Anderson et al.ll2006l : lYadavll20o4 iBellini et al] 
l2009l : iMontalto et al. I l2009l ). 



In Sec. 2 we describe the observational data and meth- 
ods used to prepare CCD images for the analysis. The pro- 
cedures employed for the determination of relative PMs of 
individual stars and for the determination of the absolute 
PM of the cluster are presented in Sec. 3. A discussion of our 
results is given in Sec. 4 and this is followed by a Summary. 



2 DATA SELECTION AND PREPARATION 

The images of M55 analyzed in this paper were selected from 
the data collected by the Cluster AgeS Experiment (CASE; 
iKaluznv et"all 120051 ') project during the period 1997-2008. 
Observations were made with the 2.5-m du Pont telescope 
at Las Campanas Observatory. All images were obtained 
with the same detector and the same set of filters. We used 
the 2K^ Tek#5 CCD camera with a scale of 0.259" /pixel 
and a field of view of 530" x 530". CASE observed M55 as 
part of a program to obtain follow-up photometric observa- 
tions of specific variables and to survey the central part of 
the cluster for detached eclipsing binaries. For the present 
analysis we selected four pointings, named F1-F4. The lo- 
cations of these fields are shown in Fig. 1. The field centers 
are at the following (RA, Dec) coordinates for Epoch 2000: 
Fl (294.9975°, -30.9621°), F2 (294.9788°, -30.9316°), F3 
(294.9668°, -30.9499°) and F4 (295.0200°, -30.9707°). The 
set of frames selected for a given observing run and for a 
given pointing is called a 'data-set'. Each data-set includes 
the V frames with the best seeing obtained for a given ob- 
serving run. We considered only frames with seeing better 
than l.l"and obtained at an air-mass less than 1.1, avoiding 
observation obtained through clouds or with a bright sky. A 
summary of the final data-sets is given in Table 1. 

For each data set we constructed an averaged, high 
signal-to-noise ratio frame using the DIAPlE package. In 
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brief, individual frames were transformed to the coordi- 
nates of a chosen reference image. The point spread func- 
tion (PSF) of each frame was also transformed to match the 
PSF of the reference image (i.e. the transformation kernel be- 
tween a given frame and the reference image was derived). 
The best seeing frame from a given data-set served as the 
reference image. The transformed frames were subsequently 
stacked by averaging to form what we later call an 'average 
frame'. During this procedure the reference as well as the in- 
dividual frames were divided into 16 overlapping sub-frames 
to reduce effects of PSF variability. 

We decided to obtain PMs for stars detected on the 
averaged frames of fields Fl and F2 for seasons 1997 and 
1999, respectively. These frames are characterized by par- 
ticularly good seeing and high S/N. We label these the 
'master' frames. The observations of fields Fl and F2 span 
9 and 11 years, respectively. Fields Fl and F2 show sig- 
nificant overlap with fields F3 and F4. Therefore, most of 
stars detected on Fl-1997 and F2-1999 images are present 
on the remaining averaged frames. A master list of stars 
was extracted from frames Fl-1997 and F 2 - 1999 using the 
DAOPHOT/ALLSTAR package fjStetsonI [l987l l. The re- 
ductions were made on 16 overlapping sub-frames result- 
ing from a 4x4 division of the respective averaged frames. 
A Mofi^at function with linear spatial variability was used 
to characterize the PSF. Due to crowding the photometry 
was extracted in an iterative way, gradually decreasing the 
detection threshold. An effort was made to avoid artificial 
splitting of bright stars which can happen when one uses an 
automatic procedure to detect missed objects in star sub- 
tracted images. During the last iteration the residual im- 
ages were inspected by eye to find previously undetected 
stars. The master lists contained 47333 and 42725 objects 
for fields Fl and F2, respectively. In the next step these 
lists were transformed to the coordinates of the averaged 
frames for different epochs. Profile photometry (as well as 
PSF modeling) was measured for these frames with ALL- 
STAR parameter set to re = 1. In this way the positions of 
objects from the master lists have been re-determined for a 
given averaged frame. As before, we used a 4x4 mosaic of 
sub-frames for actual reductions. 



3 PROPER MOTIONS 

3.1 Proper motion measurements 

Relative PMs of individual stars were derived from their po- 
sitions measured at different epochs with respect to nearby 
cluster stars. This was done using a pro cedure resembling 
that described by [Anderson et al.l (|2006l ) . As a first guess 
we selected as cluster stars objects located along the main 
sequence and red giant and horizontal branches on the CMD. 
For this purpose we constructed deep V/ {V — I) CMD's for 
both analyzed fields, including only stars with high quality 
photometry. The photometric quality was judged based on 
CHI and SHARP fit parameters returned by ALLSTAR. 
PM measurements were only made for stars with V < 21.0. 
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For any given star we selected a set of 'grid stars' lo- 
cated inside an ~80"x80"square, centered on the star to be 
measured. Only likely cluster members were considered at 
this step. The median number of grid stars was equal to 203 
for Fl and 202 for F2. The grid stars were used to calculate 
the local geometrical transformation between the appropri- 
ate master frame and a given average frame. This transfor- 
mation was calculated using the IRAF0 immatch.geomap 
and immatch.geoxytran tasks. We used a two-dimensional 
3rd order Chebyshev polynomial as the transformation func- 
tion. Subsequently, we calculated the expected coordinates 
{Xc,Yc) of the star on the average frame based on its 
(XojYo) coordinates on the master frame. The relative mo- 
tion was derived as the difference between the calculated and 
observed coordinates: dX — Xc — Xq and dY — Yc — Yq. 
The shifts dX and dY were measured for all suitable average 
frames. Finally, the PM of a given star was calculated by a 
linear fit to dX and dY as a function of time. Figure 2. shows 
as an example the PM measurement for star #13300383. 

Individual points were weighted by the amount of fiux 
corresponding to a given 'average frame'. In practice this was 
done by summing fiuxes measured for several bright isolated 
stars on individual frames before averaging. We then took 
the inverse value of the summed flux to compute the relative 
weight for each 'average frame'. PM measurements were only 
attempted for objects with positions determined on at least 
four epochs spanning at least 4 years. We retained only those 
objects for which the signiflcance of a linear fit to the PM 
exceeded 99% (i.e. a t^st was performed with a false 
detection rate of a = 1.0 %). 

PM measurements were performed in two steps. In the 
first step we made preliminary PM estimates. In the sec- 
ond step we used only reference stars for which reliable 
PMs were obtained in step one. Additionally, positions of 
reference stars were shifted back to the reference epoch (Fl- 
1997 or F2-1999) to minimize errors due to the movement 
of 'grid stars'. These new positions were used to calculate 
geometrical transformations in the second iteration. Figure 
3 shows that the total RMS of the geometrical transforma- 
tion does not change significantly with time and is weakly 
seeing-dependent . 

Equatorial coordinates for the measured stars were ob- 
tained using frame astrometric solutions b ased on a set of 
stars selected from the UCAC3 catalogue l|Zacharias et al.l 
l2010i ). We considered only stars with V < 17 which resulted 
in 671 and 636 stars for fields Fl and F2, respectively. The 
average residual of the adopted solutions is 0.15 arcsec. Our 
catalog of M55 stars includes a total of 16954 stars with a 
measured relative PM. Table 3 presents the first few lines of 
the complete catalogue which is available in electronic form 
at html : // case . camk . edu . pi. 

3.2 Error estimates 

Figure 4 shows the difference in measured PMs for 11398 
stars measured on both fields Fl and F2. For 66% of these 
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stars the difference is smaller than A/x = 0.34 mas/yr. This 
can be adopted as a robust estimate of the average error of 
a PM determination for the whole sample. For 95% of the 
stars the difference is smaller than A/i = 0.94 mas/yr. 

Figure 5 shows the dependence of logiocr^j on V magni- 
tude. For stars with V « 19.0 the median value of is equal 
to 0.47 mas/yr. At V = 20 the median value of o"^ reaches 
0.65 mas/yr. An increased scatter in ct^ can be noted for 
V < 14.5. This arises from the saturation of bright stars 
on some frames resulting in fewer available epochs for the 
PM measurement. Note that the cr^ values below V « 17.5 
arise due to systematic effects and for fainter magnitudes 
are related to photon noise statistics. The spread of the 
values may be due to the methodology of the reductions (fit- 
ting and coordinate transformations) and to the number of 
points used in the astrometric solutions. 



3.3 Completeness 

We assessed the completeness of our PM survey as a func- 
tion of V magnitude and radial distance from the cluster 
center. Completeness was defined as the ratio of the number 
of stars with final PM measurements to the number of stars 
for which we attempted to make a PM measurement. The 
results are shown in graphical form in Figs. 6 and 7. Only 
points with relative uncertainty smaller than 20% were plot- 
ted. The completeness exceeds 70% for 13 < < 17 and 
drops to 25% ai V = 20.0. As expected, the completeness 
increases with distance from the cluster center, fiattening 
at a distance of about 4 arcmin. The larger incompleteness 
observed at the cluster center is due to effects associated 
with crowding. The apparent peak at d = 4 arcmin, evident 
for all magnitude ranges, results from the fact that for more 
distant stars often only a few epochs were available. No at- 
tempt was made to estimate the completeness of the initial 
master list of stars detected in the studied field. 



3.4 Cluster membership probabilities 

A vector point diagram (VPD) including 16945 stars with re- 
liably derived relative PMs is presented in Fig. 8. Only 2.3% 
(394 objects) of the measured objects have fj, > 3.0 mas/yr 
and these are likely field stars. To estimate the probabil- 
ity that a star is a clu ster member w e use the 'local-sample 
method' described in iPlataid (Il984|}_ with further improve- 
ments proposed by iPlatais et al. In our case, how- 
ever, field stars comprise only a very small fraction of the 
total sample. 

We therefore decided to assign stars to one of three 
groups based on location in the VPD diagram and on the 
error of the measured PM. All objects with a measured PM 
were divided into five sub-samples in magnitude: 13 < V < 
17 (N = 1048), 17 < 1/ < 18 (N = 2511), 18 < V < 19 (N 
= 5140), 19 < < 20 (N = 5715) and 20 < V < 21 (N = 
2531). For each bin we selected stars with /i < 1.8 mas/yr 
and calculated for these average values and standard devi- 
ations of fj,aCosS and us {Ma, Ms, Sa, Ss)- The total value 
of the standard deviation of the PM was then calculated 
as 5 = {Sa + 5|)^^^. Average values and standard devia- 
tions of th eindividual errors of jiaCOsS and ^s were also 
calculated {MEa.s, SEa,s)- Subsequently, the total values 



ME = {MEl + ME'lf/'^ and 5*^; = {SEl + SEf^^^ were 
calculated. The results are listed in Table 2. The histograms 
showing the distribution of both components of PM for the 
brightest and for the faintest bin in V are shown in Fig. 9. 
Stars were assigned to three classes of membership. Those 
with M > 3.0 X S are considered non-members (mem=0; 
N = 379). Possible members (mem=l; N = 187) are stars 
which M ^ 3 X 5 but with > ME + 3.0 x SE). Objects 
with M < 3 X 5" and ME + S.Ox SE are considered PM 
members of the cluster (mem=2; N — 16379). The member- 
ship class defined as above is listed in 12-th column of Table 
3. 



3.5 Absolute proper motions 

Among the objects remaining on the star-subtracted aver- 
aged frames there are a number of relatively faint, unre- 
solved galaxies. These were selected by visual inspection of 
both the averaged frames and averaged frames with stars 
subtracted. The sample includes a total of 70 galaxies from 
fields Fl and F2. Eleven of these are undoubtedly relatively 
bright compact galaxies. Due to crowding of the stellar field 
we decided to extract photometry of the galaxies simultane- 
ously with photometry of the nearby stars using the GAL- 
FIT 3.0 code (Peng et aL 2010). This was done on (30x30 
or 60x60 pixel^) sub-frames centered on a given galaxy and 
extracted from the relevant averaged frame. 

For every sub-image a PSF model was constructed us- 
ing a few bright stars and IRAF's daophot routines. All but 
two of the galaxies were modeled with a 2-D Sersic func- 
tion. The fit included the closest stars, some of which were 
present in galaxy foreground. The sky background was also 
fitted. The coordinates of the galaxies derived this way were 
used to calculate their PMs relative to the cluster. This was 
done using the same procedure we applied earlier to stars. 
We excluded two galaxies from the sample since their dX 
and dY shifts showed non-linear dependence on time. We 
interpret this as due to problems with the proper deter- 
mination of their coordinates. The weighted average of the 
relative PM based on the remaining 9 galax ies is (ngCosS, 
fis) = (-3.19 ± 0.37, -9.27 ± 0.80) mas/vr. iKaluziiv et all 
(120051) in their search for cataclysmic variables in M55 iden- 
tified 7 blue stellar-like objects. One of those, M55-B1, was 
classified as QSO based on its light curve. The object is 
present in our catalogue as #23100322 with a PM value 
similar to the average PM of our sample of 9 galaxies. An- 
other blue object, #22401066 (#17550 in the lKaluznv et al.1 
2005 numbering of their reference frame), has a PM indicat- 
ing that it is a cluster member. We failed to obtain PMs for 
the re maining 5 blue stellar- like objects from lKaluznv et al.l 
(|2005t ) because they are too faint. The weighted average PM 
based on 9 galaxies and M55-B1 is {fiaCosS, /is) = (—3.30 
± 0.20, -9.14 ± 0.17) mas/yr. A VPD for these 10 objects 
is shown in Fig. 10. The point with the largest uncertainty 
corresponds to an edge-on galaxy. After dropping this object 
we obtained {^laCosS, ^is) = (-3.31 ± 0.10, -9.14 ± 0.15) 
mas/yr. This value was obtained assuming {^aCosS, fig) — 
(0.0,0.0) mas/yr. Ou r value can be compared with result of 
iDinescu et ail (|l999l ) who obtained a PM for M55 of (-1.57 
± 0.62, — 10.14± 0.64) mas/yr. Their measurement was ob- 
tained using wide-field photographic plates with a base-line 
of 25 years. 



4 K. Zloczewski, J. Kaluzny and I. B. Thompson 



4 DISCUSSION 

4.1 Color-magnitude diagram 

Figure 11 shows a.V/{B- V) CMD of M55 based on pho- 
tometry from JKIO. The left panel includes about 10 thou- 
sands stars with good quality photometry. By selecting stars 
classified as certain PM members (mem=2 in Table 3) we 
obtained the CMD shown in the right panel of Fig. 11. 

4.2 Binary stars 

A clear sequence running parallel to the main-sequence can 
be seen in the cleaned CMD. We interpret these objects 
to be binary stars in M55. The bina ry content of M55 was 
discussed at length bv lSollima et al.l (2007). These authors 
used HST/ACS images centered on the cluster core. Our 
data show that a population of binaries with g ~ 1 is also 
present in the more external parts of M55. The binary se- 
quence crosses the cluster turn-off and then extends to the 
blue including several candidate blue straggler (BS) stars. In 
addition to the BSs we note the presence of candidate yellow 
and red straggler stars. These objects are located above the 
main-sequence turnoff on both sides of the lower part of the 
subgiant branch. A definitive determination of the member- 
ship status of these red/yellow/blue stragglers can be made 
based on measurements of their radial v elocities. M5 5 has a 
large radial velocity, Vrad = 174.8 km/s (|Harrislll996l ). and it 
should be easy to separate field stars from cluster members 
with high confidence. 

4.3 Variable stars and blue stragglers 

Our PM survey can be used to assign tentative member- 
ship for variables and for other interesting objects from the 
cluster field. We co mpiled a list of var i ables and candidat e 
blue stragglers from loTech et all j 19991 ), iPvch et all l|200lh . 
iLanzoni et al] (|2007l ') & JKIO. Our PM catalog includes 52 
out of 71 known variables and 46 out of 65 candidate BSs. 
These stars are listed in Table 4. In Fig. 12 we present a 
CMD of M55 showing the upper main-sequence and BS re- 
gions. All variable BSs with available PMs turned out to 
be probable members of the cluster. The BS region con- 
tains some PM members of M55 not known to be variables. 
The light curves of these stars obtained by JKIO were re- 
examined in detail but none of them showed any evidence 
for variability. Two of the variable BSs are eclipsing bina- 
ries. Variables which with high confidence do not belong to 
the cluster are: V15, V49, V50 and V51. In addition the 
BSs B SS-7, BSS-27, BSS-31 & BSS-39 from Lanzoni et al.1 
([2003) are probable field stars. Individual cases should be 
checked taking into account individual PM errors. Further 
radial velocity studies of all variable stars are needed to re- 
liably establish their membership and evolutionary status. 

4.4 Sagittarius dSph 

The variable V15 has been found bv lOlech etal] (| 19991 ) to 
be an RR Lyr star belonging to the Sagittarius dSph galaxy. 
Part of this extended galaxy is presen t in the backgroun d 
of M55 (|lbata. Gilmore. fc Irwinlll994l : iMateo et al]fl996D . 

We measured the PM of V15 to be (naCosS, /^f) = (1.14 



± 0.31, 7.81 ± 0.36) mas/yr. The variable is located on 
the VPD inside an apparent clump formed by a few dozen 
stars. We identify this group of stars as likely members of 
the Sagittarius dwarf. For 29 stars surrounding V15 on the 
VPD we obtained a weighted mean absolute PM of {fiaCosS, 
^i4) = (-2.23±0.14, -1.83±0.24) mas/yr. This position in 
the VPD was calculated iteratively with the PM of V15 
as a starting point and using stars within a circle of ra- 
dius of r = 1 mas/yr. The measured value of the average 
PM does not change beyond the quoted errors if this pro- 
cedure is repeated for radii of 0.75 and 1.25 mas/yr. The 
positions of these stars on the cluster CMD are shown in 
Fig. 11. We conclude that this subsample is composed of 
upper main-sequence and subgiant s tars from th e Sagittar- 
ius dwarf galaxy. For a comparison, iPrvor et al.1 (|2010l ) de- 
termined the absolute PM of the Sagittarius dwarf to be 
{^iaCos6, fig) = (-2.75 ± 0.20, -1.65 ± 0.22) mas/yr based 
on archival HST images. Their measurement was made at 
positions located 7.3 & 9.7 deg NW from the center of M55. 
Earlier de terminations of the PM for the galaxy were r e- 
ported bv lDinescu et al.1 (|2005l ) and bv llbata etal] (119971 ). 



5 SUMMARY 

We have performed the first proper motion study of the 
globular cluster M55 based on ground-based CCD images. 
Relative PMs were derived for 16945 stars to V < 21.0 over 
a temporal baseline of up to 11 years. We assigned a cluster 
memberships status for all stars with available PM mea- 
surements. The absolute PM of M55 was measured to be 
{fiaCosS, fis) = (-3.31 ± 0.10, -9.14 ± 0.15) mas/yr. The 
absolute PM was derived for a portion of the Sagittarius 
dSph galaxy located in the background of the cluster. The 
results of our survey allow a selection of cluster variables, 
as well as cluster blue/yellow/red stragglers, with high cer- 
tainty. 
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Table 1. Summary of M55 observations used in this study. 
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Table 2. PM statistics for magnitude bins. Only stars with fj, < 
1.8 were used. See text for details. 
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This paper has been typeset from a TJiX/ I^T^JC file prepared 
by the author. 
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Table 3. First lines of the electronically available PM catalogue. Columns; (1) star 
ID (starting with 1 and 2 for fields Fl and F2, respectively); (2) & (3) equatorial 
coordinates (a, <5)2000.0 for epochs 1997.41 and 1999.46 for Fl and F2, respectively; 
(4) & (5) XY pixel coordinates on reference frames; (6)-(9) PMs and their errors; 
(10) number of epochs used; (11) temporal baseline; (12) cluster membership (for 
explanation see Sec. 3.4); (13) V magnitude. 
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Table 4. Tentative membership status (mem) for variables and for BSs from Lanzoni et al. (2007; IDl). 
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Figure 1. Positions of the Fl, F2, F3, F4 fields overlaid on a DSS R-band image of M55. The field of view is 15'xl5'. 
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Figure 2. An example of PM determination for one of the fast moving stars. 
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Figure 3. RMS of geometrical transformation versus time for field Fl. Average seeing for a given frame is indicated. 




Figure 4. Differences of PMs measured for stars in common to fields Fl & F2 as a function of V magnitude. 
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Figure 5. PM errors. 
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Figure 6. Completeness of the PM catalog as a function of projected distance from the cluster center. 
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Figure 7. Completeness of the PM catalog as a function of V magnitude. 
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Figure 8. The proper motion Vector Point Diagram for 16945 stars from M55. Tiie catalog includes 15 high velocity stars located outside 
of the plotted area. Stars with /i >3.0 are marked with large dots (394 objects). The arrow points to the clump of stars presumably 
associated with the Sgr dSph galaxy. 
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Figure 9. The distribution in iXcosS a^iid IJ-S for magnitude bins V = 13-17 and V = 20-21. The right panel shows a VPD for stars with 
PM < 1.8 mas/yr. Circles indicate 3xS which is defined in Sec. 3.4. 
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Figure 10. Relative PMs of 9 compact galaxies and the QSO M55-B1 measured with respect to M55 stars. 




Figure 11. Left panel - V/ {B — V) CMD of M55 based on JKIO photometry. Stars selected as likely PM members of the Sgr dSph galaxy 
are marked with squares (see Sec. 4.3).; Right panel - the CMD including only probable cluster members. Candidate blue/yellow/red 
stragglers are marked with large dots. 
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Figure 12. V/{B — V) CMD showing the turnoff and BS regions of M55. Only Ukely cluster members are plotted. Circles - pulsating 
variables of SX Phe type; squares - eclipsing binaries. 



